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Abstract: We assess the prospects for indirect detection of Weakly Interacting Massive 
Particles using a large Liquid Argon TPC detector. Signal events will consist of energetic 
electron (anti)neutrinos coming from the decay of r leptons and b quarks produced in 
WIMP annihilation in the core of the Sun. Background contamination from atmospheric 
neutrinos is expected to be low, thanks to the superb angular resolution and particle 
identification capabilities provided by the considered experimental set-up. We evaluate the 
event rates predicted for an annihilating WIMP as a function of its elastic scattering cross 
section with nucleons. This technique is particularly useful for WIMPs lighter than ~ 100 
GeV. 
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1. Introduction 

Astrophysical observations, probing gravitational potentials, provide overwhelming evi- 
dence for the existence of matter that does not emit or absorb electromagnetic radiation 
and therefore is known as dark matter (for recent reviews see Refs. @, H |). Perhaps 
the most convincing evidence for the existence of dark matter comes from the rotation 
curves of a large number of spiral galaxies Additional data, such as that collected 

by WMAP ||, favours a standard cosmological model in which the Universe is flat and 
contains about 25% non-baryonic dark matter Q. 

Notwithstanding the large set of independent observations favouring its existence, the 
exact nature of dark matter remains one of the most puzzling mysteries of our time. Par- 
ticle Physics, through extensions of the Standard Model, can provide potentially viable 
candidates for dark matter, most of which fall in the class of Weakly Interacting Massive 
Particles (WIMPs). WIMPs are stable (or very long-lived) particles, with masses ranging 
from the GeV to the TeV range, which can account for a significant fraction of the dark 
matter density in the Universe. Among the plethora of WIMP candidates, the lightest 
neutralino is the most theoretically developed |8|, ^, |l0|, pi] ]. Another possibility is offered 



by Universal Extra Dimension models [12|, in which the Lightest Kaluza-Klein particle is 
a plausible dark matter candidate [jDj ]. 

Parallel to theoretical developments, there also exists an intense and challenging exper- 



imental activity devoted to WIMP detection. So far, only a single experiment, DAMA [14], 



claims to have found evidence for the presence of WIMPs in the Galactic halo. Using a 
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direct detection method, they have measured an annual modulation, over seven annual cy- 
cles (107,731 kg day total exposure), consistent with expectations from a WIMP signature. 
Other experiments, probing similar regions of the parameter space, have found negative 
results ID], 16, 17, 18]. Therefore the situation remains highly controversial. Indirect 



detection methods have produced so far negative results |T^, p0[ . 

In this paper we concentrate on indirect detection methods. We assess, in a model- 
independent way, the capabilities that a massive Liquid Argon (LAr) Time Projection 
Chamber (TPC) offers for identifying neutrino signatures coming from the products of 
WIMP annihilations in the core of the Sun. Unlike existing measurements and future 



searches in large Cerenkov neutrino telescopes [21, 22, 24], where a statistically signif- 
icant excess of high energy muons is expected to be measured, here we propose to look 
for an excess of electron-like events over the expected background of cosmic neutrinos. 
In particular, our signature will be a high energy v e and v e charged current interactions 
pointing in the direction of the Sun. We evaluate our sensitivity and present the expected 
event rates over the mass region of 10 to 100 GeV where this technique is most effective. 
Note that this includes the mass range favoured by the DAMA experiment [14|. 



2. Neutrinos from WIMP Annihilation in the Sun 

WIMPs that constitute the halo of the Milky Way can occasionally interact with massive 
objects, such as stars or planets. When they scatter off of such an object, they can po- 
tentially lose enough energy that they become gravitationally bound and eventually will 
settle in the center of the celestial body. In particular, WIMPs can be captured by and 
accumulate in the core of the Sun. 

The capture of heavy particles in the Sun was first calculated in Ref. [25]. Improved 



formulae for WIMP capture were later given in Ref. [^6| . WIMPs can accrete as a result of 
both spin-dependent and spin-independent (scalar) interactions. The solar WIMP capture 
rate is given by |2t| : 

C® ~ 3 SSxlO 2 ^- 1 ( Plocal ^ ( — km / s \ 3 f g H,SD + gH,si + 0-07otte,sA ZlOOGeVV 
X S V0.3 GeV/cm 3 / V ^locai / V 10" 6 pb J \ m wm p J 

(2.1) 

where pi OC ai is the local dark matter density, v\ oca \ is the local rms velocity of the halo 
WIMPs, (Th,sd> <7h,si and <JH e ,si are the spin-dependent and spin- independent WIMP- 
Hydrogen and WIMP-Helium cross sections, respectively, and mwiMP the mass of the 
WIMP. Throughout this paper, we will refer to the quantity cth.sd + o"h,si + 0.07o"He,si 
simply as the elastic scattering cross section, cr e lastic- 

In this paper, we do not consider the possibility of observing neutrinos from WIMPs 
accumulated in the Earth. Given the smaller mass of the Earth and the fact that only 
scalar interactions contribute, the capture rates for our planet are not enough to produce, 
in our experimental set-up, a statistically significant signal. For a detailed discussion of 
WIMP capture in the Earth, see Ref. [§7| • 

Accumulated WIMPs will annihilate into ordinary matter (quarks, leptons, gauge and 
Higgs bosons). In the case of neutralino dark matter, annihilations do not directly produce 
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neutrinos, as they are Majorana particles || . Instead, neutrinos with a spectrum of energies 
below the mass of the WIMP will be generated in the decays of annihilation products. 
Other dark matter candidates may produce neutrinos directly, however. For example, 



Kaluza-Klein dark matter [28, 29] or sneutrino dark matter |3(], |3lJ can each annihilate 
directly to neutrino pairs. 

As already pointed out, we are only interested in the expected fluxes of electron neutri- 
nos and anti-neutrinos. To compute this, we use the analytic expressions given in Ref. p2| . 
The differential energy flux is: 

where Ta is the WIMP annihilation rate in the Sun and R is the distance from the Sun to 
the Earth. The sum is over annihilation channels, Bp being the fraction of annihilations 
which go to channel F. (dN/dE)p is the differential energy spectrum of electron neutrinos 
at the surface of the Sun per WIMP annihilating through channel F. 

The final electron (anti) neutrino flux is obtained using the following set of assumptions: 

• We take typical values for the local dark matter density and local rms velocity, namely 
Piocal = 0.3 GeV/cm 3 and Vi oca i = 270 km/s. 

• Equilibrium has been reached and therefore the annihilation rate, Ta, is equal to the 
capture rate in the Sun, C®. This is a reasonable assumption given the range of cross 
sections we consider here. 

• We restrict ourselves to a low mass WIMP scenario, namely to-wimp < 100 GeV; in 
this way we minimize neutrino flux depletion due to propagation effects across the 
Sun. 

• Due to the former assumption, phase space considerations avoid WIMP annihilation 
into top quarks and Higgs bosons. Although in the narrow mass range between 80 
and 100 GeV, annihilations to gauge boson pairs Z°Z°) are possible, we 
do not explore this channel. We instead focus on WIMP annihilations to the bb and 
t + t~ channels. These are generally the dominant channels for a neutralino in this 
mass range. 

The electron neutrino spectrum from annihilations to r + r~ is given (in the relativistic 



limit) by |3j 



dN\ = / 2r r _^ p \ _ 3x2 + 2x3) e _ E/Ek ^ 

dE ) T+T _ V mwiMP / 



where x = E/mwiMP and E is the neutrino energy. r r — *evv — 0.18, the branching 
ratio of taus into an electron plus neutrinos. The factor, e~ E ^ Ek , accounts for the 
absorption of neutrinos as they propagate through the Sun. E^ ~ 130 and 200 GeV 



for u e and v e , respectively [33 
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The electron neutrino spectrum from annihilations to bb is similar 



m 



(—) = ( )(1 - 3x 2 + 2x 3 ) e- E ' E * (2.4) 
\dE J bl V^^wimp/ 

Here, x = E/(zm-wiMp) where z = 0.7 is the fraction of energy retained after 
hadronization. The inclusive branching ratio for b quarks to electron neutrinos, 
Fb->evX, is approximately 0.11. 

• We consider that tau pairs are produced with a branching fraction B T + T - = 10% 
while for b quarks we assume a branching ratio of B b i = 90%. 



3. The Liquid Argon TPC 

The detection device we propose to use for indirect WIMP detection is a Liquid Argon 
Time Projection Chamber (LAr TPC). This technology was first put forward in 1977 p5| . 
In a LAr TPC, the ionization charge, released at the passage of charged particles, can 
be transported over large distances thanks to the presence of a uniform electric field. The 
signals induced by drift electrons are recorded by a set of successive anode wires planes thus 
providing three dimensional image reconstruction and calorimetric measurement of ionizing 
events. Unlike traditional bubble chambers, limited by a short window of sensitivity after 
expansion, the LAr TPC detector remains fully and continuously sensitive, self-triggerable 
and without read-out dead time. 

The feasibility of the LAr TPC technology has been demonstrated thanks to an ex- 
tensive R&D programme, developed by the ICARUS collaboration, that culminated in the 
construction and operation of a 600 ton detector |3(|. A surface test of this detector has 
demonstrated the high level of maturity reached by this technology |37], |38|, |4(], f|l]]. 

Future applications of this detector technology have been proposed to study neutrino 
properties within the context of neutrino factories ]^2|, |43| or to study neutrinos from 
supernovae [H. An interesting long-term application is the possibility to build a giant 
LAr TPC detector with a mass of 100 ktons. The technical aspects and physics prospects 



of such detector have been extensively discussed elsewhere [45, 46]. In the present paper 
we evaluate for the first time the capabilities that such a detector will have in the search 
for dark matter WIMPs. 

We have focused our interest on v e , v e charged current events, since a LAr TPC offers 
very good electron identification capabilities. In addition, all final state particles can be 
measured and identified (see figures |l| and [2] where a full simulation of a v e charged current 
interaction is shown). Thanks to the high granularity and superb imaging capabilities, 
we expect to have a very accurate measurement of both energy and incoming direction. 
Therefore our signal will correspond to a neutrino interaction that points in the direction 
of the Sun and contains an energetic primary electron. Possible background sources are 
extensively discussed in the next section. 
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Figure 1: Simulated v e interaction in a Liquid Argon detector (longitudinal view). The electron 
neutrino [E v — 26.3 GeV) travels from left to right and suffers a charge current interaction with 
the atoms of the medium, resulting in a highly energetic electron (23 GeV), two pions (Ekin = 0.7 
and 1.4 GeV) and a proton (Ekin = 1.6 GeV). This last particle is not visible in this projection. 
The neutrino direction is computed combining the information coming from the electromagnetic 
shower and the hadronic jet. 




Figure 2: Rotated view of the same v e interaction shown in figure |l|. In this projection the four 
particles coming out from the main interaction vertex are clearly visible: two pions, one proton and 
the electron (that appears as a big black spot). 



4. Sources of v e Background 

We have identified several sources that could, a priori, significantly contribute to the ex- 
pected background rate. The most obvious is atmospheric neutrinos, where we expect 
contributions from both v e jv e charged currents and vjv neutral currents (where a 7T° is 
misidentified as an electron). In the energy range we study here (~10 to 100 GeV) solar 
neutrinos do not significantly contribute to the background. We have also considered the 
flux of neutrinos originating from cosmic ray interactions with the interstellar medium in 
the galaxy. As shown in Ref. [47], for the range of energies considered in this study, this 



flux is negligible compared to the atmospheric flux. The diffuse flux of v e from active 
galactic nuclei is also small. 

Fluxes of v T are also a concern given that r leptons, produced in charged current 
interactions, can subsequently decay into electrons. According to Ref. |4^1 , the expected 
flux of v T from sources like atmospheric neutrino oscillations, the galactic plane and galaxy 
clusters are small and will not give a significant contribution in our energy range. 

Therefore, the most important source of background is by far due to atmospheric neu- 
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trinos. The evaluation of this flux <j) v has been performed using the calculations published 
in Ref. [E9]. In this paper the authors present a detailed 3-dimensional calculation of the 
atmospheric neutrino flux based on the FLUKA Monte Carlo model |)0| . The final FLUKA ta- 
bles for three experimental sites (Kamiokande, Gran Sasso and Soudan) are also accessible 



in Ref. 151 



Figure ^ shows the differential atmospheric electron neutrino flux as a function of the 
neutrino zenith angle and the neutrino energy expected at the latitude of Gran Sasso []5l| . 
Similar distributions are built and used for D e , and v^. Only neutrinos in the energy 
range 1 - 100 GeV are shown. 

Without loss of generality, the spectra at Gran Sasso has been taken as reference for 
all calculations presented in this work. The conclusions of the paper remain valid for 
other underground laboratories such as Kamiokande or Soudan since the neutrino fluxes 
at these locations are expected to be similar (at high energies, differences are less than a 
percent f5l|). 




, 0.8 

2 04 



Figure 3: Differential atmospheric electron neutrino flux as a function of the neutrino energy and 
zenith angle. The data, taken from Refs. 49 |5l[], correspond to the Gran Sasso latitude and is 
taken as our reference to compute the expected background rate. 
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5. Data Analysis 



The number of expected background events contributing from each neutrino species (N Ui ) 
is determined from simple integration of the associated atmospheric differential neutrino 
flux (4>Ui) taking into account the involved cross sections (cr Ui ), the detector mass (M^) and 
the data taking period: 



where the integration is done above a minimum energy threshold, E™ m , and inside a 
limited solid angle on the sky, AQ, in order to reduce the background contribution. In the 
following we assume a detector mass of 100 ktons of LAr and 10 years of data taking. 

We have taken into account not only Charged Current (CC) events for v e and v e 
interactions but also Neutral Current (NC) events for both electron and muon species. 
Neutral current events contribute to the background mainly through (1) electrons from 
Dalitz decays and (2) misidentified pions. In this work we have assumed a conservative 
90% efficiency on the rejection of NC events, even though Monte Carlo studies on LAr 
detectors show much better capabilities (>99%). and CC interactions are excluded 
from the calculation since the presence of a muon in the final state would immediately 
veto the event as a possible v e (the muon-electron misidentification probability in a LAr 
detector is almost zero). 

In general, when looking at neutrinos coming from astrophysical point-like sources, 
the atmospheric neutrino background can be reduced to a very low level by selecting small 
angular regions of the sky. This applies to our particular problem where we are interested 
in neutrinos coming only from the Sun. Therefore the region of interest is of the order of 
the angular size of the Sun and, in principle, the neutrino background coming from this 
narrow sky window will be small. 

At this stage we have to estimate Af2, the apparent size of the Sun "seen" by the 
detector. Its value will be driven by our capabilities to determine the incoming neutrino di- 
rection, i.e. the detector angular resolution. This calculation has been carried out through 
a detailed and realistic simulation of neutrino interactions inside a LAr detector. Neutri- 
nos of different species in a wide energy range (< 100 GeV) and the final state particles, 
produced after the v interaction, were tracked using the FLUKA package. As an example, 
figures [| and |2| show a picture of one of these simulated events where a 26.3 GeV electron 
neutrino (flying from left to right) undergoes a CC interaction giving rise to several sec- 
ondary particles: v e n — > e~ p ir + tt~ (2n3j). The three charged hadrons appear as dark 
thin lines over light grey background, the electron develops as an electromagnetic shower 
and the neutral particles are invisible. The reconstruction of the v direction is done using 
the information coming from all particles produced in the final state. 

The result of the simulation is shown in figure ||, where the difference between the sim- 
ulated and reconstructed neutrino directions, A9, is plotted as a function of the incoming 
neutrino energy (solid curve). Note that for the considered energy range, O(10 GeV), the 
smearing introduced by the Fermi motion of the initial state nucleon and by re-interaction 
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of hadrons inside the nucleus are negligible and the improvement on the neutrino direction 
measurement achieved by combining the electron and the hadronic jet informations is sig- 
nificant. Due to the excellent reconstruction capabilities offered by the LAr technique, the 
angular resolution remains below 2° for the highest energies up to about 15 GeV, increasing 
by a factor of 4 at lower energies (~ 5 GeV). 

In case of an ideal detector with perfect angular resolution, the sky integration window 
would be limited to the size of the Sun (~ 6.7 • 10~ 5 sr). In our more realistic scenario Af2 
has been estimated by smearing this number with the detector resolution. The dashed line 
in figure |] shows the resulting half size of the Sun (radius) after the smearing (see right 
axis). For comparison, the real Sun radius seen from the Earth (~ 0.256°) is also plotted. 

Concerning the signal, the number of expected events is computed in a model inde- 
pendent way as described in section [2] by integrating equation 2J1 between E™ tn and the 
WIMP mass. Contributions from electron neutrino oscillations in the range of considered 
neutrino energies (above 10 GeV) are negligible. Therefore, no loss or gain on signal events 
from oscillated neutrinos is taken into account. 

The final estimation for a 50 GeV WIMP mass is presented in figure [5] where we show 
the number of expected signal (squares) and background (dots) events as a function of the 
cut on the minimum neutrino energy. The numbers are normalized to a 100 kton LAr 
detector and 10 years of data taking and for a WIMP elastic scattering cross section in 
the Sun of 10 -4 pb (the result can be easily rescaled to other exposures or cross sections). 
As expected, the background rate quickly decreases when increasing the neutrino energy 
cut because of (1) the nature of the atmospheric neutrino flux and (2) the better detector 
angular resolution. From this figure we conclude that the search can be safely considered 
as "background free" for values of E™ m above ~ 10 GeV (contamination ~ 0.2 events). 

Figure || gives the evolution of the signal events as a function of the WIMP elastic 
scattering cross section, for neutrinos of energy above 10 GeV and three values of the 
WIMP mass. The background (horizontal line) is below the signal for values of Ceiastic 
above ~ 5 • 10~ 6 pb. 



6. Prospects for Discovery 

In section |E| we have evaluated the average number of signal and background events (called 
hereafter /is and hb, respectively) expected in the detector for a certain exposure, £ = 
Mass x time = 100 ktons x 10 years. At this point we would like to know whether a 
numbers of expected events from a particular dark matter model will allow the experiment 
to claim a discovery at a certain significance level. This problem is fully discussed in 



Ref. [52], where a general and well defined criterion of discovery is proposed for statistical 



studies of this nature: 

- Usually a 5a criterion (called "method A" hereafter) is used, so that in order to 
consider that a discovery will take place, the signal must be above a 5a fluctuation 
of the background: fis > 5J]Tb. 
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Figure 4: Detector angular resolution (A0) as function of the incident neutrino energy obtained 
from a Monte Carlo simulation in a LAr detector (solid line and left axis). The figure also shows 
the size of the Sun (radius) for an ideal detector (horizontal solid line) and the resulting apparent 
size after detector smearing (dashed line, read on the right axis). 



This commonly used rule frequently gives misleading (if not clearly wrong) results. 
Only when /is and \xb are large enough to describe with Gaussian distributions can 
this rule be applied. This is clearly not the case in our study. 

- Instead, we have followed the method proposed in Ref. j52j ("method B" hereafter), 
where the discovery criterion is given by two numbers, namely: e, the significance 
required to consider an excess of events a discovery, and 5, the probability that 
the experiment will obtain such an excess. For instance, when 5 = 0.5 and e = 
5.733 x 10~ 7 , we can state that at least 50% of the time, the experiment will observe 
a number of events which is 5cr or more above the background. 

The results based on "method B" (5 = 0.5 and e = 5.733 x 10 -7 ) are presented in 
figure [7[ The plot shows the minimum number of years required to observe a WIMP 
positive discovery signal as function of the WIMP elastic scattering cross section. The 
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Figure 5: Expected number of signal and background events in a 100 kton LAr detector after 
10 years of data taking as function of the cut on the minimum neutrino energy. Signal events 
correspond to a 50 GeV WIMP mass and are normalized to a elastic scattering production cross 
section in the Sun of 10 -4 pb. The number of atmospheric neutrino events includes contributions 
from muon and electron (anti)neutrinos, and takes into account the detector angular resolution. 



numbers are given for three values of the WIMP mass and only neutrinos above 10 GeV 
are considered. 

From this result we conclude that, in 10 years of data taking, a 100 kton LAr detector 
would be able to claim a clear discovery signal according to "method B", provided <7 e i as ti c 
is above ~ 10 -4 pb for WIMP masses greater than 20 GeV. On the other hand, only one 
year of data taking would be enough to observe a discovery signal provided <r e i as ti c is above 
5 x 10~ 4 pb and mwiMP ~20 GeV. 

Finally, we have studied the sensitivity of the results with the cut on the minimum 
neutrino energy. We remind the reader that the value of E" im sets the energy above which 
all events are integrated. By lowering the cut we would accept too many atmospheric 
background events and by increasing it both, fj,$ and \ib decrease, but at different rates 
(see figure ||) . The quantitative result in terms of discovery prospects are shown in figure || 
for two values of E™ m and a reference WIMP mass of 50 GeV. 
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Figure 6: Expected number of signal and background events as a function of the WIMP elastic 
scattering production cross section in the Sun. The three lines correspond to three values of the 
WIMP mass and a cut on the minimum neutrino energy of 10 GeV. 



7. Comparison to Other Indirect Detection Techniques 

Using neutrino detectors to search for evidence of dark matter has been studied in the 
context of several experiments. Currently, high energy neutrino telescopes that use large 
volumes of a natural Cerenkov medium (such as ice or water) show the greatest promise 
for dark matter detection. The AMANDA experiment, located at the South Pole, has 
been taking data for several years. Its successor, IceCube is currently under construction 
at the same site. When completed, IceCube, will constitute the largest neutrino detector 
to date, with a full cubic kilometer of instrumented volume. The ANTARES experiment 
is currently under construction in the Mediterranean. Although its effective volume will 
be considerably smaller than IceCube's, ANTARES has been designed to be sensitive to 
neutrinos of lower energy, and therefore may be competitive with IceCube for detecting 
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Figure 7: Discovery potential: Minimum number of years required to claim a discovery WIMP 
signal from the Sun in a 100 kton LAr detector as function of elastic for three values of the WIMP 
mass. See text for details on the applied statistical criteria. 



dark matter. 

At the energies we are considering here (well below a few TeV), neutrino telescopes 
observe neutrinos through energetic muons produced in charged current interactions. The 
direction of these neutrinos can only be reconstructed to ~ 1.2° / (TeV) , the angle at 
which the muon and neutrino are aligned. The angular cuts which can be imposed when 
searching for a signal from the Sun are thus considerably weaker than for the liquid Argon 
technique described in this article. Overcoming the atmospheric neutrino background is 
an essential requirement if a liquid Argon TPC detector is to succeed. 

Experiments such as AMANDA, ANTARES and IceCube are (or will be) considerably 
larger than the 100 kton liquid Argon TPC detector we use here. The target mass of 
IceCube is approximately 10 4 times greater, for example. For very massive WIMPs (heavier 
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Figure 8: Detector discovery potential for two values of the minimum neutrino energy cut and a 
50 GeV WIMP mass. See text for details on the applied statistical criteria. 



than 100-200 GeV), high-energy neutrino telescopes provide a considerably more sensitive 
probe of dark matter annihilations than the technique described here. For lighter WIMPs, 
however, this may not be the case. 

Considering a 50 GeV WIMP, for example, we expect neutrinos to be generated with 
energies on the order of 15-20 GeV and lower (assuming the annihilation modes discussed in 
section §). Although such neutrinos are well within the energy range that a liquid Argon 
detector could be sensitive to, they would generate muons with energies of 8-10 GeV, 
which are below the threshold for planned high-energy neutrino telescopes. Neutrinos 
from WIMPs lighter than ~ 100 GeV are likely to be missed by high-energy experiments, 
but are ideal for a liquid Argon TPC detector. 

For discussions of the prospects for the detection of dark matter with high-energy 
neutrino telescopes, see Refs. [[5^, 54, 55, 56]. 
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8. Conclusions 



In this paper we have studied the ability of a giant Liquid Argon Time Projection Chamber 
(TPC) to detect dark matter indirectly via observation of neutrinos produced in WIMP 
annihilation in the core of the Sun. We have adopted an approach which differs from 
previous studies, where an excess of upward- going muons would signal the presence of 
WIMPs in the core of the Sun or the Earth. Our method takes advantage of the excellent 
angular reconstruction and superb electron identification capabilities Liquid Argon offers 
to look for an excess of energetic electron (anti) neutrinos pointing in the direction of the 
Sun. The expected signal and background event rates have been evaluated, in a model 
independent way, as a function of the WIMP's elastic scatter cross section for a range of 
masses up to 100 GeV. 

We have also presented a prospective study that quantifies the detector discovery 
potential, i.e. the number of years needed to claim a WIMP signal has been discovered. 
With the assumed set-up and thanks to the low background environment offered by the 
LAr TPC, a clear WIMP signal would be detected provided the elastic scattering cross 
section in the Sun is above ~ 10" 4 pb. 

Acknowledgments 

This work has been supported by the CICYT Grant FPA2002-01835. S.N. acknowledges 
support from the Ramon y Cajal Program. D.H. is supported by the Leverhulme Trust. 
TW. is supported by DOE grant DE-FG05-85ER40226. A.B. and T.W. thank the Be- 
nasque Center for Physics where the early stages of this work started. 

References 

[1] C. Munoz, Int. J. Mod. Phys. A 19, 3093 (2004) [arXiv:hep-ph/0309346]. 

[2] G. Bcrtonc, D. Hooper and J. Silk, arXiv:hcp-ph/0404175. 

[3] T. J. Sumner, Living Rev. Relativity 5, 4 (2002); 

http://relativity.livingreviews.org/Articles/lrr-2002-4/download/index.html 

[4] M. Persic, P. Salucci and F. Stcl, Mon. Not. Roy. Astron. Soc. 281, 27 (1996) 
[arXiv:astro-ph/9506004]. 

[5] Y. Some and V. Rubin, arXiv:astro-ph/0010594. 

[6] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) 
[arXiv:astro-ph/0302209]. 

[7] W. L. Freedman and M. S. Turner, Rev. Mod. Phys. 75, 1433 (2003) 
[arXiv:astro-ph/0308418] . 

[8] H. Goldberg, Phys. Rev. Lett. 50, 1419 (1983). 

[9] L. M. Krauss, Nucl. Phys. B 227, 556 (1983). 



- 14 - 



[10] J. R. Ellis, J. S. Hagelin, D. V. Nanopoulos, K. A. Olive and M. Srednicki, Nucl. Phys. B 
238, 453 (1984). 

[11] P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schclke and E. A. Baltz, JCAP 0407, 008 
(2004) [arXiv:astro-ph/0406204]. 

[12] T. Appclquist, H. C. Cheng and B. A. Dobrescu, Phys. Rev. D 64, 035002 (2001) 
[arXiv:hep-ph/0012100]. 

[13] G. Servant and T. M. P. Tait, Nucl. Phys. B 650, 391 (2003) [arXiv:hcp-ph/0206071]. 

[14] R. Bcrnabci et al. [DAMA Collaboration], Phys. Lett. B 480, 23 (2000); Riv. N. Cim. 26 
(2003) 1 [arXiv:astro-ph/0307403]. 

[15] D. S. Akerib et al. [CDMS Collaboration], arXiv:astro-ph/0405033; Phys. Rev. D 68, 
082002 (2003) [arXiv:hcp-ex/0306001]. 

[16] G. Chardin et al. [EDELWEISS Collaboration], Nucl. Instrum. Mcth. A 520, 101 (2004). 

[17] N.J.T. Smith et al., talk given at the 4th International Workshop IDM2002, York, England; 
http: / /www. shcf.ac.uk/physics/idm2002/talks/pdfs/smith_n.pdf 

[18] B. Ahmed et al, Astropart. Phys. 19, 691 (2003) [arXiv:hcp-ex/0301039]. 

[19] M. Ambrosio et al. [MACRO Collaboration], Phys. Rev. D 60, 082002 (1999) 
[arXiv:hep-ex/9812020] . 

[20] S. Desai et al. [Super-Kamiokande Collaboration], arXiv:hep-ex/0404025. 

[21] J. Ahrcns et al. [The IceCube Collaboration], Nucl. Phys. Proc. Suppl. 118, 388 (2003) 
[arXiv:astro-ph/0209556]; A. Goldschmidt [IceCube Collaboration], Nucl. Phys. Proc. 
Suppl. 110 (2002) 516. 

[22] F. Halzcn et al. [AMANDA Collaboration], Phys. Rcpt. 307, 243 (1998) 
[arXiv:hep-ex/9804007] . 

[23] F. Blanc et al. [ANTARES Collaboration], Presented by L. Thompson on behalf of the 

ANTARES Collaboration at the 28th International Cosmic Ray Conferences (ICRC 2003), 
Tsukuba, Japan, 31 Jul - 7 Aug 2003. Published in Tsukuba 2003, Cosmic Ray 1743-1746. 

[24] S. Bottai [NESTOR Collaboration], Prepared for 2nd International Workshop on the 
Identification of Dark Matter (IBM 98), Buxton, England, 7-11 Sep 1998. 

[25] W. H. Press and D. N. Spergel, Astrophys. J. 296, 679 (1985). 

[26] A. Gould, Astrophys. J. 321, 571 (1987); Astrophys. J. 388, 338 (1992). 

[27] J. Lundberg and J. Edsjo, Phys. Rev. D 69, 123505 (2004) [arXiv:astro-ph/0401113]. 

[28] H. C. Cheng, J. L. Feng and K. T. Matchcv, Phys. Rev. Lett. 89, 211301 (2002) 
[arXiv:hep-ph/0207125]. 

[29] D. Hooper and G. D. Kribs, Phys. Rev. D 67, 055003 (2003) [arXiv:hep-ph/0208261]. 

[30] T. Falk, K. A. Olive and M. Srednicki, Phys. Lett. B 339, 248 (1994) 
[arXiv:hep-ph/9409270]. 

[31] D. R. Smith and N. Weincr, Phys. Rev. D 64, 043502 (2001) [arXiv:hcp-ph/0101138]. 

[32] G. Jungman and M. Kamionkowski, Phys. Rev. D 51, 328 (1995) [arXiv:hep-ph/9407351]. 



- 15 - 



[33] I. F. M. Albuquerque, L. Hui and E. W. Kolb, Phys. Rev. D 64, 083504 (2001) 
[arXiv:hep-ph/0009017]. 

P. Crotty, Phys. Rev. D 66, 063504 (2002) [arXiv:hep-ph/0205116]. 

C. Rubbia, "The Liquid-Argon Time projection Chamber: a new concept for Neutrino 
Detector", CERN-EP/77-08, (1977). 

S. Amcrio et al. [ICARUS Collaboration], Nucl. Instrum. Mcth. A 527, 329 (2004). 

S. Amoruso et al. [ICARUS Collaboration], Nucl. Instrum. Mcth. A 523, 275 (2004). 

S. Amoruso et al, Nucl. Instrum. Meth. A 516, 68 (2004). 

M. Antonello et al, Nucl. Instrum. Mcth. A 516, 348 (2004). 

S. Amoruso et al. [ICARUS Collaboration], Eur. Phys. J. C 33, 233 (2004) 
[arXiv:hep-ex/0311040]. 

Nucl. Instrum. Meth. A 508, 287 (2003) [Erratum-ibid. A 516, 610 (2004)]. 

A. Bucno, M. Campanclli and A. Rubbia, Nucl. Phys. B 589, 577 (2000) 
[arXiv:hcp-ph/0005007]. 

A. Bucno, M. Campanelli, S. Navas-Concha and A. Rubbia, Nucl. Phys. B 631, 239 (2002) 
[arXiv:hcp-ph/01 12297]. 

I. Gil-Botclla and A. Rubbia, JCAP 0310, 009 (2003) [arXiv:hcp-ph/0307244]; JCAP 0408, 
001 (2004) [arXiv:hcp-ph/0404151]. 

A. Rubbia, Proceedings of the Second NO- VE International Workshop on Neutrino 
Oscillations in Venice; pages 321-350. Venice, Italy, 2003. Ed. M. Baldo-Ceolin. 

A. Ereditato and A. Rubbia, arXiv:hep-ph/0409143. 

G. Ingelman and M. Thunman, arXiv:hep-ph/9604286. 

H. Athar and C. S. Kim, arXiv:hep-ph/0407182. 

G. Battistoni, A. Ferrari, T. Montaruli and P. Sala, Astr. Part. Phys. 19, 269 (2003) 
[arXiv:hep-ph/0207035]. 

See http : //www . f luka . org/ for further details. 

FLUKA flux tables are available at http : //www.mi . inf n. it/^battist/neutrino .html 

J. J. Hernandez, S. Navas and P. Rebecchi, Nucl. Instrum. Mcth. A 372, 293 (1996). 

L. Bergstrom, J. Edsjo and P. Gondolo, Phys. Rev. D 58, 103519 (1998) 
[arXiv:hcp-ph/9806293]. 

V. D. Bargcr, F. Halzen, D. Hooper and C. Kao, Phys. Rev. D 65, 075022 (2002) 
[arXiv:hep-ph/0105182]. 

J. L. Feng, K. T. Matchcv and F. Wilczek, Phys. Rev. D 63, 045024 (2001) 
[arXiv:astro-ph/0008115]. 

V. Bertin, E. Nczri and J. Orloff, Eur. Phys. J. C 26, 111 (2002) [arXiv:hcp-ph/0204135]. 



- 16 - 



